Gluino Stransverse Mass 



Won Sang Cho 1 , Kiwoon Choi 1 , Yeong Gyun Kim 1 ' 2 and Chan Beom Park 1 
1 Department of Physics, KAIST, Daemon 305-017, Korea 
2 ARCSEC, Sejong University, Seoul 143-747, Korea 

We introduce a new observable, 'gluino stransverse mass', which is an application of the Cambridge 
niT2 variable to the process where gluinos are pair produced in proton-proton collision and each 
gluino subsequently decays into two quarks and one LSP, i.e. gg — > qqxi IQXi- We show that the 
gluino stransverse mass can be utilized to measure the gluino mass and the lightest neutralino mass 
separately, and also the 1st and 2nd generation squark masses if squarks are lighter than gluino, 
thereby providing a good first look at the pattern of sparticle masses experimentally. 
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PACS numbers: 14.80.Ly, 12.60.Jv 

The Large Hadron Collider (LHC) at CERN will soon 
explore TeV energy scale, where new physics beyond the 
Standard Model (SM) likely reveals itself P, . Among 
various new physics proposals, weak scale supersymmc- 
try (SUSY) [3(] is perhaps the most promising candidate, 
which provides a solution for gauge hierarchy problem 
and complies with gauge coupling unification. Further- 
more, with R-parity conservation, the lightest supersym- 
metric particle (LSP) becomes a natural candidate for 
the non-baryonic dark matter in the Universe. 

Once SUSY signals are discovered through event excess 
beyond the SM backgrounds in inclusive search chan- 
nels, the next step will be the measurements of SUSY 
particle masses and their physical properties in various 
exclusive decay chains. Then it might be possible to 
reconstruct SUSY theory, in particular SUSY breaking 
mechanism from the experimental information on SUSY 
particle masses. In this regard, determination of gaug- 
ino masses has a significant implication because predic- 
tions of the gaugino masses are rather robust compared 
to those on sfermion masses 

In this paper, we introduce a new observable, 'gluino 
stransverse mass', which is an application of vari- 
able [B| to the process where gluinos are pair produced 
in proton-proton collision and each gluino subsequently 
decays into two quarks and one LSP, 



VP 



gg -> qqx° qqx°, 



(1) 



where q stands for the 1st or 2nd generation quark. We 
show that the gluino stransverse mass can be utilized to 
measure the gluino mass and the lightest neutralino mass 
separately, and also the 1st and 2nd generation squark 
masses if squarks are lighter than gluino, thereby provid- 
ing a good first look at the pattern of sparticle masses 
exp er iment ally. 

If light enough, gluinos would be pair produced co- 
piously in proton-proton collision (pp — > gg), and each 
gluino decays into two quarks and one LSP (g — > qqXx) 
through three-body decay induced by an exchange of off- 
shell squark or two body cascade decay with intermedi- 
ate on-shell squark. For each gluino decay g — > qqxi, a 



transverse mass is constructed, which is defined as 



z / qq x 

i T (m qq x,m x , pf,p T 



l qqT T III. 
P « • p* j 



(2) 



where m qq T and are the transverse invariant mass 
and transverse momentum of the qq system, respectively, 
while m x and p T are the assumed mass and transverse 
momentum of the LSP, respectively. The transverse en- 
ergies of the qq system and of the LSP are defined by 



P§?l a 



Px\ 2 



(3) 



With two such gluino decays in each event, the gluino 
stransverse mass 711x2(9) is defined as 



i T2 



(9) 



max{? 



2(!) 2( 2 h 



} , (4) 



where the minimization is performed over all possible 
splitting of the observed missing transverse energy lss 
into two assumed transverse momenta, p£^ and p^ 2 ^. 

From the definition ((4]) of the gluino stransverse mass, 
one obtains a relation: 



m T2{g) < n^g for TO} 



(•5) 



i.e. 171x2(9) is l ess than or equal to the true gluino mass 
m g if the trial LSP mass m x is equal to the true LSP 
Therefore one can determine ra g from the 
endpoint measurement of 111x2(9) distribution: 



m^(m x ) 



max [m T2 (g)] , 
all events 



(G) 



if we already know the true LSP mass. However m 
might not be known in advance and then, for a trial m x 
which is different from m^o, m™2 X ( m x) wm differ from 
Trig and can be considered as a function of m x . 

As we will see, m^^^x) nas different functional form 
depending upon whether the 1st and 2nd generation 
squarks are heavier or lighter than gluino. Thus, in order 
to investigate the TO x -dependence of to™| x , we consider 
two cases separately. 
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If squark masses are larger than gluino mass, the gluino 
will undergo three body decay into two quarks and an 
invisible LSP through the exchange of an off-shell squark: 



9 



qqxi- 



(7) 



In order to get an idea on how my| x (m x ) are determined 
for a generic value of the trial LSP mass m x , we consider 
two extreme momentum configurations and then con- 
struct the gluino stransverse mass associated with each 
of them. 

The first momentum configuration is that two gluinos 
are produced at rest and then each gluino subsequently 
decays into two quarks moving in the same direction, and 
one LSP whose direction is opposite to the quark direc- 
tion. Furthermore, two sets of gluino decay products are 
parallel to each other and all of them are on the trans- 
verse plane with respect to the proton beam direction. 
For such momentum configuration, we have 



l qqT 



™% = 0, 



(8) 



where the final state quarks are regarded to be massless. 
The transverse energies and transverse momenta of the 
qq systems are given by 



7.99(1) 



7-99(2) 
J T 



2m g 



E^(max), 



\pf 2) \ 



(9) 



where m g and to x o are the true gluino mass and the 
true LSP mass, respectively, and the corresponding total 
missing transverse momentum is 

\p™ iss \ = 2E T q (max). (10) 

It has been shown that for certain momentum config- 
urations, the rriT2 variable ('balanced solution') can be 
obtained as the minimum of subject to the following 
two constraints Q: 



Tflrp j 



Pt + Pt 



X(2) 



(11) 



If one applies this balanced solution approach to the mo- 
mentum configurations with the above ro^y , E T q ^ 1,2 \ 
p9?(i,2)^ pmiss anc j g^jj undetermined p^ 1 ' , one finds 
that the minimum of rrSp{= nij') is obtained when 



p^ 1 ) = p* { z > = p™ lss /2, leading to the following gluino 
stransverse mass: 



m T2 {g) = E T q (max) + yj ' (E T q (max)) 2 + m| (12) 

for generic m x . One can show that this 111x2(9) corre- 
sponds to to™| x for m x < m x o [y]: 



m^(m x ) = 



to? 



2mg 



\ 



to- — min 
9 x" 



2ms 



for m x < to-^o. 



(13) 



Note that to™2 X ( to x ~ m x?) = m 9 as required- 
Other extreme momentum configuration which would 



m^l x (m^ 



for to x > m x o is that gluinos are 



determine 

pair produced at rest and for each gluino decay, two 
quarks are back to back to each other while LSP is at 
rest. In addition, all particles are on the transverse plane. 
In this case, one easily finds 



m 



(2) 
qqT 



TOo 



m q q T (max), (14) 



and also 



with 



Ef {l) = E q T q(2) = m qqT (max), 



Pt 



rmss 

Pt 



(15) 



(16) 



For the momentum configurations with m ' T , 



E T q(ia \ p™ (1 ' 2) and p T liss given by dHJ), (JH]), (16]) re- 
spectively, m!j) is equal to m^ for all possible split- 
ting of the missing transverse momentum: p™ ss = = 

p T ^ + p T ^ 2 \ and the minimum of to^'(= ml ) occurs 
p^ 2 '' = 0. Then the gluino stransverse 
mass obtained as a balanced solution is given by 



when p T 



mT2(g) = m qq T(max) + m x 



(17) 



for generic value of m x . This in fact corresponds to m™f J 



for m x > m x o 



111 



T2 



: (to x ) = [m- g -m x o) +: 



for 



l x > TO x o 



(18) 



which again gives to™| x (to x 



a) 



The above momentum configurations leading to the 
expression (fT8|) have p™ ss = 0, thus could be eliminated 
by the event cut imposing a lower bound on |p™ ISS | in 
the real data analysis. However, a more detailed study 
shows that there exist momentum configurations yielding 
the same expression of for m x > m Xi , while 

having a sizable \p T llss \ comparable to tos/2, for instance 
a configuration in which the two quarks from the first 
gluino move in the same transverse direction, while the 
other two quarks from the second gluino are back to back. 
As a result, the functional behavior of (|18|) can still be 
constructed from collider data under a proper cut on the 
missing transverse momentum. 

By now, it should be clear that m™* for m x < m x o 
(Eq. (fT3|) ) has a quite different form from to™2 X for 



> m x o (Eq. JIB])). As required, they should cross 
at m x = to x o. Thus, if the function to™ 2 1x (to x; ) could be 
constructed from experimental data, which would iden- 
tify the crossing point, one will be able to determine the 
gluino mass and the LSP mass simultaneously. 
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The experimental feasibility of measuring rrig and m^o 
through TO™2 X depends on the systematic uncertainty as- 
sociated with the jet resolution since m™f x is obtained 
mostly from the momentum configurations in which some 
(or all) quarks move in the same direction. Our Monte 
Carlo study indicates that the resulting error is not so sig- 
nificant, so that rrig and rn^o can be determined rather 
accurately by the crossing behavior of m™| x . As a spe- 
cific example, we have examined a parameter point in the 
minimal anomaly mediated SUSY-breaking (mAMSB) 
scenario [7] with heavy squarks, which gives 

m g = 780.3 GeV, m x « = 97.9 GeV, 

and a few TeV masses for sfermions. We have gener- 
ated a Monte Carlo sample of SUSY events for proton- 
proton collision at 14 TeV by PYTHIA 8]. The event 
sample corresponds to 300 integrated luminos- 

ity. We have also generated SM backgrounds such as 
it, W/Z + jet, WW/WZ/ZZ and QCD events, with less 
equivalent luminosity. The generated events have been 
further processed with a modified version of fast detec- 
tor simulation program PGS which approximates an 
ATLAS or CMS-like detector with reasonable efficiencies 
and fake rates. 

The following event selection cuts are applied to have 
a clean signal sample for gluino stransverse mass: 

1. At least 4 jets with P T i,2.3,4 > 200,150,100,50 
GeV. 

2. Missing transverse energy E™ lss > 250 GeV. 

3. Transverse sphericity St > 0.25. 

4. No b-jets and no leptons. 

For each event, the four leading jets are used to calcu- 
late the gluino stransverse mass. The four jets are di- 
vided into two groups of dijets as follows. The highest 
momentum jet and the other jet which has the largest 
\pj e t\AR with respect to the leading jet are chosen as 
the two 'seed' jets for division. Here pj e t is the jet mo- 
mentum and AR = y A<j> 2 + Arj 2 , i.e. a separation in 
azimuthal angle and pseudorapidity plane. Each of the 
remaining two jets is associated to a seed jet which makes 
the smallest opening angle. Then, each group of the di- 
jets is considered to be originating from the same mother 
particle (gluino). 

Fig[T] shows the resulting distribution of the gluino 
stransverse mass for the trial LSP mass m x = 90 GeV. 
The blue histogram corresponds to the SM background. 
Fitting with a linear function with a linear background, 
we get the endpoint 778.0 ±2.3 GeV. The measured edge 
values of 1x1x2(9), i- e - to t2 X j as a function of m x is shown 
in FigO Blue and red lines denote the theoretical curves 
of (fl3|) and IT8|) . respectively, which have been obtained 
in this paper from the consideration of extreme momen- 
tum configurations. (A rigorous derivation of (| 1 3[) and 



(TTSf will be provided in the forthcoming paper @.) Fit- 
ting the data points with the curves (fT3|) and (fl"8|) , we 
obtain m g = 776.3 ± 1.3 GeV and m x o = 97.3 ± 1.7 GeV, 
which arc quite close to the true values, m g = 780.3 GeV 
and m x o = 97.9 GeV. This demonstrates that the gluino 
stransverse mass can be very useful for measuring the 
gluino and the LSP masses experimentally 
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FIG. 1: The rriT2(g) distribution with m x = 90 GeV for the 
benchmark point of mAMSB with heavy squarks. Blue his- 
togram is the SM background. 
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FIG. 2: m™!* as a function of the trial LSP mass m x for the 
benchmark point of mAMSB with heavy squarks. 

Let us now consider the case that squark mass rrig is 
smaller than the gluino mass rrig. In such case, the fol- 
lowing cascade decay is open; 

g -> qq^ qqx°i- (19) 

In this case also, we consider two extreme momentum 
configurations which are similar to those considered for 
three body gluino decay, and construct the corresponding 
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gluino stransverse masses. Here again, we assume gluinos 
are pair produced at rest and each gluino decays into a 
quark and a squark on the transverse plane with respect 
to the proton beam direction. 

If the quark from squark decay is produced in the same 
direction as the first quark from gluino decay, and the two 
sets of gluino decay products are parallel to each other, 
the transverse energies and transverse momenta of the qq 
systems are given by 



H/rp Hirp 



2m g 



= |p«W| = |p«W| 



Ej?(max), 



(i) 

qqT 



m q 2 qT f° r the transverse invariant 



and we have m 
masses of the qq systems and |p™ ss | = 2E T q (max) for 
the total observed missing transverse momentum. Then 
the same procedure to obtain the gluino stransverse mass 
(fT2| can be applied to this case, leading to m™2 X which 
is same as Eq. flU]) : 



m^(m x ) 



m~ 

Xi 



2Ulg 



\ 




for m x < rn^o. 



(20) 



Now we consider other extreme momentum configura- 
tion in which the quark from squark decay is produced in 
the opposite direction to the first quark from gluino de- 
cay. In this case, the invariant transverse masses, trans- 
verse energies and transverse momenta of two qq systems 
are given by 



b qqT 



99(1) 



e: 



Ip t 9(1) 



2(2) 
l qqT 

99(2) 



(m?-m?)(m?-m|o) 



E, 



<?9(2) 

Pt 



m 
m 



-(1 



-(1 



-(1 



o 

xi 



m 



with the total missing transverse momentum given by 

remiss 

Pt 



-2p» ,(1) . 



Imposing the two constraints, m T 



(i) 



rrvp and 



~miss 

Pt 



p T ^ + p T ^ 2 \ on the momentum configura- 



tions having mffl , E T q(1 ' 2) , pf h2) ~ miss 



, p ™^(= -2p™ (1) ) 



as above, we obtain the following balanced solution of 

„x(2) _ _„99(1). 

Pt Pt 



m T i{g) at p T (1) 



m T2 {g) 



2(1) . 

qqT + 

2{E q T q(1) 



(21) 



p«(l)|2 + m 2 + | p «(l)| 2) 



Again, one can show that this represents rn™^ for m x > 



m x1 Pi' yielding 



i T2 



m 
m 



9 

ml, 

x 



'-) 



\ 



m 5 



-(1 



(22) 



Note that the two functions (|2"0"|) and (|2"2"|) cross at to x = 
m^o for which m™| x = m^. 

If one could construct (|2"D|) and (|22p accurately enough 
from data, one would be able to determine all involved 
sparticle masses, i.e. nig, m^o and uiq. To see how feasi- 
ble it is, we examined a parameter point of mirage media- 
tion model [1Q], providing the following sparticle masses; 



m 5 



,, - 821.4 GeV, 



694.0 GeV, m^o = 344.2 GeV. 



We have generated a Monte Carlo sample for this 
benchmark point of mirage mediation, corresponding to 
100/6 -1 integrated luminosity. After the event selection 
cuts similar to the case of three body gluino decay, we 
obtain FigJ3] showing the distribution of my 2 (ff) for the 
trial LSP mass m x — 350 GeV. The edge value as a 

function of m x is shown in Fig0J Fitting the data points 
to the curves ([20]) and ([22]), we obtain m g = 799.5 ± 11.1 
GeV, mg = 678.2 ± 7.0 GeV and m^o = 316.7 ± 15.4 
GeV. Though the overall scale of the fitted values are 
well close to the true sparticle masses, the central values 
are somewhat lower than the true values, which is mainly 
due to a mild crossing of two curves. The situation can 
be improved if we include the information from squark 
stransverse mass for the process pp — > qq — > qXiQXii 
providing a relation between the edge value of the squark 
stransverse mass and the trial LSP mass 61: 



m™| x (squark) 



m? 
x" 



2niq 



\ 



m1 
x" 



2m?, 



Including such information, we get m g — 803.4±6.0 GeV 
and m^o = 322.4 ± 7.7 GeV for the benchmark point. 
The discrepancy between the fitted mass values and the 
true mass values may still come from various systematic 
uncertainties such as the effects of event selection cuts, 
which is beyond the scope of this paper. 

To conclude, we have introduced the gluino stransverse 
mass, and shown that it can be used to determine the 
gluino mass and the lightest neutralino mass separately, 
and also the 1st and 2nd generation squark masses if 
squarks are lighter than gluino. 
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FIG. 3: The m,T2(g) distribution with m x = 350 GeV for the 
benchmark point of mirage mediation. 













Li 


ght squark 















































o ion zoo mu 40<) son 600 



FIG. 4: m^f* as a function of the trial LSP mass m x for the 
benchmark point of mirage mediation. 
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